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upper bound for the deepest edge state. In the case of Abelian gauge theories we analyse 
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of the bounds on the boundary conditions and the size of the system is derived under 
general assumptions. 
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1. Introduction 

Quantum fields are the fundamental pillars of high energy physics. In the last few 
years, it has also been shown that many new condensed matter systems have an 
effective description in terms of relativistic quantum fields. This novel perspective 
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not only was useful for developing new experiments, but also to understand some 
open questions in field theory and many body physics. In particular, the lack of 
analytic tools to study strongly correlated quantum systems and the difficulties 
of some numerical approaches can now be overcome by using condensed matter 
analogue systems. Moreover, they can be used as quantum simulators to understand 
some non-perturbative effects which were elusive by standard techniques. 

One of the main differences between condensed matter and fundamental particle 
physics is that in condensed matter, the materials where the effective field theory 
lives are of finite size and have boundaries, in contrast with the unbounded nature 
of Minkowski space-time. This difference is crucial and brings some extra features 
to the former like the appearance of a particular type of bound states (edge states) 
which are localised at the boundary of the system. These states give rise to new 
effects and transport phenomena. In particular, they are responsible for the appear¬ 
ance of edge currents and edge conductivity in graphene and topological insulators 
[28]. Physical effects where edge states also play a relevant role include the quantum 
Hall and Casimir effects. 

A specially interesting case is when edge states appear at the boundary between 
two phases, e.g. metal-superconductor or metal-insulator. The contact between a 
gapless phase and gapped one induce interesting dynamical effects on the boundary 
which are enhanced by the quantum theory. 

In this paper we address the characterisation of edge states in field theories 
defined on Riemannian manifolds M with metric g and regular, codimension one, 
boundaries dM. The metric gives volume forms dpg on M and d/ig on dM and their 
associated Hilbert spaces "Hm and 'Hom ■ In the case of bosonic field theories (scalar 
or gauge theories), the operator which accounts for the modes associated with the 
kinetic energy comes from the Laplace-Beltrami operator —Am- Consistency of 
the quantum field theories requires that this operator is self-adjoint on Hm and 
positive, otherwise the quantum vacuum is unstable and unitarity of the quantum 
theory is lost. See [6,7,22] and references therein in this connection. In the Euclidean 
space M = , these properties are guaranteed. The operator —Am is essentially 

self-adjoint and positive on smooth functions of compact support in 

However, on manifolds with boundaries, it does not always have the required 
properties. In particular, for scalar fields in the domain Dm of the Laplacian, 

(d>i, —AM‘f’2) = / d^g AMd’2) 

Jm 

= [ dAigWl-V$2- [ dMg$l|9^V„$2|aM ’ (1-1) 

JM JdM 

where, as mentioned above, dM is the induced manifold at the boundary, dpg and 
d/ig denote the volume forms on M and dM respectively and n is the outward- 
drawn normal on dM. Notice that the second term, unlike the first, can be negative 
for $1 = $ 2 . 

In the case that the boundary is an interface boundary between two phases, we 
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just consider it as a double boundary: one face for each side with opposite normal 
derivatives, and again asymmetric boundary conditions on the two faces can give 
rise to negative contributions. 

It has been known for a long time [10,37] that such a negative spectrum does 
occur for the Robin boundary conditions 

V„$ = Md>, (1-2) 

for positive values (/i > 0) of the Robin parameter and also for more general bound¬ 
ary conditions [26]. In particular, Asorey et al. [8] proved that, when /r becomes 
large and the Dirichlet condition is approached, the negative spectrum recedes to 
—oo and at the same time the corresponding eigenstates get progressively more 
localised at the boundary. The existence of this phenomenon is shown in [7,29] for 
even more general boundary conditions. Numerical procedures for analysing the 
associated eigenvalue problem are also described therein. 

The physical interest for these boundary conditions relies in the fact that they 
arise in a natural way at the interface between materials in different phases [4]. 
The best known examples come from electromagnetism between the normal and 
superconductor phases. The non-interacting order parameter field $, that in general 
can be a tensor field, is determined by the spectrum and eigenvectors of —A -|- 
instead of those of —A. The shift in the spectrum provided by the term m? is such 
as to lift the negative energy levels to non-negative values. At the same time, the 
bulk levels, which already have positive eigenvalues for the case m = 0, acquire a 
gap of order m^. Topological insulators also share these features. 

From a mathematical viewpoint the situation above requires the relevant oper¬ 
ator —A to be lower bounded. This is to ensure that the shift by a mass term 
is able to make the operator —A -|- positive. Lower boundedness of self-adjoint 
extensions of elliptic operators on compact manifolds was proved by G. Grubb [26]. 
The case of non-compact domains with compact boundaries was proved very re¬ 
cently by the same author [27]. Lower boundedness of general self-adjoint extensions 
of non-elliptic operators like d*d on one-forms, which describes the electromagnetic 
field, remains an open problem. 

This article is focused on the study of the general structure of edge sates for 
field theories with scalar, vector and spinor fields. Edge states are associated with 
the negative energy eigenstates of the relevant operators, cf. [7,29,4]. 

First we present an alternative proof of the well-known semiboundedness of the 
Laplace-Beltrami operator in a compact Riemannian manifold with boundary for 
Robin boundary conditions, cf. [26,19,20]. This is done in Section 2. The proof 
is based on a recent approach [32] that can be used for more general boundary 
conditions than those considered here and that works also in the context of non¬ 
compact manifolds with compact boundaries. As pointed out there, generalisation 
to non-compact boundaries is also possible. The advantage of this proof is that 
it keeps track of the different geometrical structures that naturally appear in the 
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problem and can be used to provide also upper bounds for the ground states. 
This is done in Section 2.2. The approach provides estimations of the shape of 
the eigenfunctions while the upper bounds coincide with the asymptotic results 
obtained in [42]. When they exist, the large negative energy states are located at 
the boundary earning the name of edge states. 

Furthermore, the proof can be modified to treat the situation of the non-elliptic 
operator d*d on one-forms mentioned above. The literature has no result on this 
problem despite its physical relevance for the electromagnetic field. This is done in 
Section 3. In order to prove lower boundedness of this non-elliptic operator, one has 
to take into account the infinite dimensional kernel of the operator. This is done 
by means of a convenient use of gauge invariance. The bounds do not depend on 
the gauge chosen and are thus completely general. The situation for fermionic fields 
is considered in Section 4 and Section 5. We analyse the edge effects arising from 
different choices of the boundary conditions. In this case the relevant operator is 
the Dirac operator whose spectrum is unbounded below and above. Edge states are 
associated in this case with the appearance of eigenstates with eigenvalue A within 
the mass gap, i.e. —m < X < m. Again, a convenient modification of the proof 
in Section 2 is used to show the existence of such states for appropriate boundary 
conditions. In Section 4 we consider Atiyah-Patodi-Singer boundary conditions. 
Surprisingly, in contrast with what happens for scalar fields and vector fields, there 
is a threshold size of the manifold below which the edge states disappear. This 
leads to qualitatively different behaviour in this situation. In Section 5 chiral bag 
boundary conditions are considered showing similar results. 

2. Edge States in Scalar Theories 

Let us consider for simplicity a complex massive scalar free field (j) with mass m on 
a manifold M with a compact regular boundary dM. 

The quantum dynamics of the field is affected by the presence of the boundary. 
The main effect is due to the contribution to the kinetic term of the quantum 
Hamiltonian governed by the quadratic differential operator — A -|- m^. Boundary 
effects appear through the boundary conditions that test fields have to satisfy to 
preserve unitarity of the quantum theory [7]. The explicit requirement is that — A-|- 
rn^ has to be a self-adjoint, positive operator. This consistency condition imposes 
severe constraints on the type of boundary conditions that stable physical systems 
must satisfy. 

We shall consider the family of general Robin boundary conditions 

V„$(p) = p{p)^{p) , p G dM , (2.1) 

where p is any smooth and hence bounded function which is not necessarily positive. 
This includes the superconducting case, where /r = m > 0, as a special case (see 
[4]). Although these are not the most general local boundary conditions that one 
can consider, see [30], they are general enough to present the phenomena of edge 
states associated to the presence of boundaries. 
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Among the eigenstates <I>„ of — A satisfying this boundary condition for positive 
/i(p), there are states with negative eigenvalues, i.e., — A<I)„ = Xn^n with A„ < 0. 
As ii{p) approaches infinity, the boundary condition approaches Dirichlet boundary 
condition = 0, the negative eigenvalues approach minus infinity and the 

corresponding eigenvectors get localised near dM (see [8] for details). 

However, the theory may still be consistent and the vacuum stable because of 
the presence of the mass term. Even if —A has a negative spectrum, the whole 
physical operator —A + may be positive. This requires the spectrum of —A 

induced by the boundary condition function /i(p) to be bounded from below by 
2 

—m. 

Although, mathematically speaking, p{p) might be any non-trivial function on 
the boundary dM, when modelling the behaviour of a superconducting phase, the 
fact that the penetration length in the superconductor is small across the whole 
boundary {skin effect) requires that /i(p) > 0 for any p € dM. Notice that in any 
case the boundary condition (1.2) is invariant under time-reversal as long as m and 
fj,{p) are real. Hence, like topological insulators, the system is T-invariant, with edge 
excitations and an incompressible bulk [4]. 


2.1. Lower bound theorem for Robin boundary conditions 

From now on and throughout the rest of the article there will appear several different 
Hilbert spaces associated to differentiable manifolds. We will denote the space of 
square integrable functions over a Riemannian manifold S as (S). The associated 
scalar product and norm will be given respectively by (•, •)e: and || • Hh . Incidentally, 
the subindices may be dropped when referring to the manifold M. In addition, we 
will need to use the Sobolev spaces of order I and 2 that we shall denote as 'H^(S) 
and 'H^(S) respectively. These spaces are the natural spaces to define differential 
operators of first and second order. Their corresponding norms will be denoted as 
II • ||«i(B) and II • ||?^ 2 (B) respectively. We refer to [3] for further details on these 
spaces. 

We shall show now that there exists a lower bound, —/Xg, for the spectrum of the 
Laplace-Beltrami operator with Robin boundary condition defined on a Riemannian 
manifold M with Riemannian metric g. The result is given by the following theorem, 
which implies that the quantum field theory is consistent whenever > /Xg . 

Theorem 2.1. Let M be an oriented Riemannian manifold with regular oriented 
boundary dM. Let n be the outgoing normal vector and /x a smooth function, 
both defined at the boundary dM . The Laplace-Beltrami operator —A restricted 
to smooth functions satisfying the Robin boundary conditions 

V„$(p) = p{p)^p) , pGdM , (2.2) 

is essentially self-adjoint and bounded from below, i.e. 


(2.3) 
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with Ho a finite positive constant. 

Although the existence of the lower bound (2.3) is known [26], for completeness 
we include the proof here. Our alternative proof is based on cobordism methods 
which will be very useful to obtain upper bounds in the next section and can be 
extended to non-elliptic operators in the case of gauge theories in Section 3. 

In the proof we shall make an intensive use of the quadratic forms associated to 
the operators, cf. [44,18]. While there is a one-to-one correspondence between semi- 
bounded self-adjoint extensions and quadratic forms, the quadratic forms have the 
advantage that while maintaining the semibounds they can be defined on domains 
larger than those of the operators themselves. More importantly, in some cases even 
the boundary conditions are suppressed by passing to the quadratic form. This hap¬ 
pens, for instance, in the well-known case of Neumann boundary conditions where 
the domain of the associated quadratic form is the full Sobolev space of order 1, 
without restrictions, cf. [18, Theorem 7.2.1]. 

Let us finally remark that for negative Robin parameters the associated self- 
adjoint extensions of the Laplace-Beltrami operator are trivially positive. Indeed, 
the right hand side of (1.1) is positive for $1 = $2 . On the contrary, for positive 
Robin parameters there is no straightforward way to prove semiboundedness. The 
reason is that the £^(9M)-norm of the trace, i.e. of the restriction to the bound¬ 
ary of a given function at the bulk, cannot be bounded by the £^(M)-norm of the 
function. This is one of the conclusions of the well known Lions Trace Theorem, cf. 
[39]. The proof of the necessary and sufficient conditions for semiboundness of even 
order elliptic operators appearing in [26,27] relies essentially on the fact that the 
inverse of the operator with Dirichlet®' boundary conditions is a compact operator. 
This clearly does not hold for the non-elliptic operator mentioned above and thus 
this proof cannot be used in this case. 

Let us consider Riemann normal coordinates in a collar neighbourhood S of the 
boundary dM. Let r be the radial normal coordinate increasing from r = Rq — e 
at the inner boundary of the collar to r = Rq at the physical boundary dM of the 
system. In these coordinates, 

E=[Ro- e, i?o] X dM, (2.4) 

while the Riemannian metric g of M restricted to S can be written as 

(o5(r,0)) ’ 

in terms of boundary coordinates 6 := {0*} of dM. 

Due to the splitting of M into S and its complement M\S in M (see Fig. 1), 
a new boundary appears. We denote it by (9S_ . Notice that the boundary of the 
collar is therefore dE = dM U dE-. This boundary has to be considered whenever 


Strictly speaking it is the so called hard extension of the operator. 
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Fig. 1. Decomposition of the manifold M into the collar neighbourhood E and its complement. 
0E_ defines a new boundary. 

one defines differential operators on S and M\S. We consider different Hilbert 
spaces on each manifold associated to this splitting and denote them as follows. 
The spaces of square integrable functions are going to be denoted as 
and £^(M\2). The Sobolev spaces of order 1 will be: ^^(5) and . 

These are going to be the natural spaces to define the quadratic forms. Similarly, 
the Sobolev spaces of order 2 will be: and . The scalar 

products and norms associated to these spaces are determined by the restrictions 
(pull-backs) of the Riemannian metric g and distinguished by subscripts. We will 
also need to consider three different Laplace-Beltrami operators —Am, —Ah and 
—Am\h each one defined on a different domain I?m, and • 

Proof of Theorem 2.1. That the boundary conditions (2.2) lead to an essentially 
self-adjoint operator is a well-known result, see for instance [25,8,32]. The only 
restrictions concern the regularity of the function , but we consider that it is 
smooth. 

The proof of the semiboundedness can be derived in 5 steps. 

i) The collar Laplacian. Let — Ah be the Laplace-Beltrami operator on the collar 
5 with metric defined by the restriction of g. Consider as its domain the space 
of functions that satisfy the Robin boundary condition at dM and the Neumann 
boundary condition at the auxiliary boundary , i.e. 

Ph = {$ e ^2(5) I v„$(p) = pip)^ip) ,pedM-, v„$(p) = 0 ,p G a5_}. 

( 2 . 6 ) 

The operator defined this way is essentially self-adjoint, 
ii) Bound on the collar Laplacian. 

The quadratic form associated to the collar Laplacian of the previous step 
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and therefore —are bounded from below, i.e., 

($,-Ab$)b >-/ioll^lll. e . (2.7) 

The constant po is of the order of \pip) \ ■ The bar denotes the closure 

with respect to the corresponding norm. We postpone the proof of this step until 
the end of this section. 

iii) The inner Laplacian. 

Consider the Laplacian —Am\e with Neumann boundary condition on , 
i.e. defined on the domain 

t ^ m\e = {$ e I v„$(p) = 0 e 9 S_} . (2.8) 

This operator is also essentially self-adjoint. It is clearly non-negative since 
Green’s identity gives 

($, -Am\h^)m\s = (V$ , V$)m\h > 0 

Moreover, [18, Theorem 7.2.1] ensures that the quadratic form associated to this 
operator has domain 'H^(M\S) without any constraints. 

iv) Equivalence of the bounds of the operator and of the associated 
quadratic form. 

We want to find a lower bound for the operator — Am defined on 
Pm = {$ e n^M) I Vn$(p) = p{pMp) ,p€dM}, 
that is, to find a non-negative constant po such that 

($,-Am$) > 

The relation between self-adjoint operators and the associated quadratic forms 
establishes that this is equivalent to showing that 

($, -Am$) > ^ e , 

where ||| • |||am stands for the graph-norm of the quadratic form defined by 

P|||i^:= (cf>,AM<i>)+m||d>f , 

where m is a lower bound of the quadratic form. 

A sufficient condition for this inequality to hold is 

($, -Am$) > ^ e . 

This is so because the graph-norm of the quadratic form is continuous with respect 
to the Sobolev norm of order one, i.e. 

|||‘h|||AM A Ar||<l>||^i(M) ■ 

To prove the latter inequality it is enough to perform integration by parts once 
and use the boundary condition to get 


($ , —Am^) = (V<I> , Vd)) — (<I)|aM , p{p)^\dM}dM 

<liv$f+ sup|Mp)II1$UII^m- 
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The first term is clearly bounded by the norm of the Sobolev space of order 1. 
The second one is bounded by means of the Lions trace inequality, cf. [3,39], 

||‘b|aM||aM < • 

v) Sum of quadratic forms. 

First notice that because of [18, Theorem 7.2.1], the Neumann boundary con¬ 
dition at the auxiliary boundary disappears when taking the closures with 
respect to the Sobolev norm of order 1 and therefore we have that 

^ll-lkyM) ^^II-IIkPh) . 

The set on the left hand side is strictly contained in the set of the right hand side 
because functions on the right hand side are allowed to be discontinuous along 
the auxiliary boundary while this is not the case on the left hand side. Now 
for any $ € , we have that 

($, = ($ , -Ah$)h + ($ , -Am\h^’)m\h (2.9a) 

(2.9b) 

where the quadratic forms in (2.9a) are those defined in Steps i) and iii). The 
inequalities of (2.9b) follow by Step ii) and the fact that 


Proof of Step ii). 

Given a fixed parameter 0 < <5 < 1, the parameter e defining the width of the 
collar neighbourhood S can be chosen such that the bounds 

{l-S)^\g{Ro,e)\ < y/WM < il + S)^\~g{Ro,e)\ (2.10) 


hold. 

Now let I := [i?o — £,Ro] ■ According to the boundary conditions of Eq. (2.6), 
we have 


($ , -Ah<I))2 = 


/ / ^ (d^>, d^)dgg_ - {g} , 

JdM JI 


r /• 5 $ 

laMJi dr dr 


dHg_ + g ^(d6i$,de$)d^g_ - {g},ii(p)oM 

JdM Jl 


> 




fdM JI 

>{l-6) 


IdM 


MI 


5<j) ,9$ 

dr dr 


dr 


l-,5 




dPg 


( 2 . 11 ) 


where ip = ^\gM denotes the boundary value of $, while g- := gls, R = sup/r(p) 
and de<l> denotes the components of d$ that are tangent to the boundary. Notice 
that there is only a boundary term contribution on dM because the Neumann 
boundary condition cancels the boundary contribution of the inner boundary . 
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The first inequality holds because the second term of the second line is positive 
definite. 

The last inequality above ensures that a lower bound of (<&, — can be 

given by a lower bound of the quadratic form on the right hand side. It is easy to 
show that this latter quadratic form is the quadratic form associated to the operator 

d^ 

-A^:=-—(g)I (2.12) 

dr^ 

defined on C^{I) ® C^{dM) by the collar radial projection of the Laplacian. The 
operator is densely defined on C^{I) and is essentially self-adjoint with mixed 
boundary conditions: Neumann boundary condition at {i?o ~ e} and Robin bound¬ 
ary condition 

at i?o • Moreover the lower bound of — A^ coincides with that of — gyj- For a more 
detailed proof we refer to [32], where more general boundary conditions are consid¬ 
ered. 

In order to obtain the bound (2.7), it is therefore enough to show that the 
one-dimensional problem of the Laplacian with Robin boundary condition at one 
boundary and Neumann boundary condition at the other is bounded from below. 

Consider the following one-dimensional spectral boundary value problem for 
functions </> in the interval [0, ej: 

-(j)" = E<j), ()i'(0) = 0 <i)\e) = c<l){e), OO. (2.13) 

We can obtain its spectrum by finding the zeros of the spectral function that is 
obtained from the plane wave ansatz ^(x) = -I-The spectral func¬ 

tion for any type of boundary condition in [0, e] has been analysed in [6,7]. In the 
present case, it reduces to 

fctan(efc) = —c . (2-14) 

Therefore there is an infinite number of positive eigenvalues = k^. However, 
there is a special bound state with negative energy Eq = — which is obtained 
from the spectral equation (2.14) assuming that k = iK is imaginary and k is given 
by the solution of the equation 

«:tanh(eK) = c . (2-15) 

This equation has only one solution for positive c. Thus we have shown that 
(2.12) is semibounded. Therefore the quadratic form on the right hand side of (2.11) 
is semibounded with the same constant. This shows (2.7) and therefore Theorem 2.1. 

□ 

The relations above show the dependence on the various parameters of the 
problem. 
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Fig. 2. Bounds on the negative energy edge mode for a particle on [0,1] with Neumann boundary 
condition at 0 and Robin boundary condition ^^(1) = c0(l), c > 0 . The red curve is the exact 
solution and the upper and lower bounds correspond to Km(c) = + tanh c ) + tanh c ) 

and km(,c) = ctanhc. 


The solution k(c) of (2.15) can be found numerically [29], but it can also be 
bounded from above, k(c) < Km(c) j by the following function: 


Km{c) 


(fh + 


ec 

tanh €c 


) tanh(^ + 


—'^1 
tanh ec ' 


(2.16) 


(See Fig. 2 where we set e = 1). 

This bound provides an explicit analytic bound for the quadratic form: 


($, -Ah$)h > 





(2.17) 


i.e., Ho = Km(i^)- 

Let us try to give some further estimates. Assume that en ^ 1. Then 
Ktanh(e«:) « k and thus the spectral condition becomes k ~ c whenever eK « 
ec 1. In the case that c = h/{1 — 5), this gives 


Ho = K 


H 

{1-S) 


whenever ck ^ 1, that is, whenever h /{1 — 5) « /io ^ b 

Now assume that eK ^ 1. In this case we have that «:tanh(eH;) « ck^ and thus 
the solution of the spectral equation becomes k « -y/c/e whenever en « <C 1. 
Thus, we have that 


Ho 


H 


e{l-S) 


H 

l-(5 



whenever 
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2.2. Existence of edge states with Robin boundary conditions 

In the previous section we have shown that for Robin boundary conditions the 
Laplace-Beltrami operator is always bounded from below. We will show now that 
when the Robin parameter is positive, p > 0, negative energy edge states appear 
for large positive values of p. This result has been proved in [8] in a more general 
framework so that it applies to our particular case. However the bounds that we 
are going to obtain are sharper. To prove the existence of such a negative energy 
eigenstate it will be enough to obtain a negative upper bound for the expectation 
value function of the Laplace-Beltrami operator for some function in an appropriate 
space. 

From Green’s identity it is straightforward to obtain the upper bound 

($ , — Am$)m = (v$ , V$)m — (v? , P‘p}dM < (v$ , V$)m — p{^ , ‘p}dM, (2.18) 

with $ € Vm and where p := inf p{p). This implies that the quadratic form associ¬ 
ated to the Laplace-Beltrami operator — Am with local Robin boundary condition 
given by p{6) is bounded above by the form associated with the Robin boundary 
condition with constant parameter p. 

We shall need the following slight generalisation of Davies’ Theorem, [18, Theo¬ 
rem 7.2.1]. It states that smooth functions with any value of the normal derivatives 
at the boundary are dense in , cf. [32, Lemma 4.1]. 

Lemma 2.1. Let 4) € and f G Then, for every e > 0 there exists 

$ G C°°{M) such that ||$-$|li < e, [jv? -'^ll«i/2(aM) < e and ||/-<^||«i(om) < £• 

As in the proof of Theorem 2.1, it is enough to obtain the bound (2.18) for 
some function in $ G T>m^^ Hwh") ignima above establishes that T>m^^ _ 

TL^{M). We can use again the splitting of the manifold showed in Fig. 1 and for 
functions $ G 77^(2) 0 77^(M\S) we can split the quadratic form as 


($ , -A„$)m = ($ , -A^$)h 0 ($, -A^^^^)m\s ■ (2.19) 

Now, consider functions on (() G Ve which only depend on the radial coordinate 
r, 4>{r), and extend them to the whole M by defining 


,x = (r,0)G 
[())(i?o-e) ,xGM\S 


We have that 

($H , -Am$h)m = , -Ah$h)h 

JdMJRo-e 9r dr 
<{l + d) 


— p{4>S , 4>B.)dM 


f 

O 

1_ 

'd^E d^E' 

JdM 

J Ro-e 

-dr dr - 


( 2 . 20 ) 


dpg. ( 2 . 21 ) 
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Thus, since the quadratic forms at both sides are dehned on the same domain, 
it is enough to hnd an upper bound of the one-dimensional problem (2.13) 
of the Laplacian on the interval [0, e] with Robin boundary condition at r = 0 and 
Neumann boundary condition at r = e. 

As we have seen in the previous section, this one-dimensional operator has an 
edge state $0 which is given by the zero of the spectral function (2.15) with Robin 
parameter c = /i/(l -I- 5). The solution of equation (2.15), k(c) , is bounded from 
below, k.m{c) < k(c), by 

«^m(c) := ctanhec. (2.22) 

Thus if we take d)H(r, 6) = <i>o(r) for (r, 0) G S we get: 

, -A„<I>h)m < - [«M(Tf3)] Vh , <I>h)h, (2.23) 

which shows the existence of negative modes of the Laplace-Beltrami operator. 
Notice that even though the function km(c) is positive, the eigenvalue associated 
to the boundary condition defined by c is —km{c)^ ■ These modes are edge states 
and are localised at the boundary of M. Notice that this upper bound for edge states 
is lower than the bound found in Ref. [8], although the latter one applies to more 
general boundary conditions. Every parameter, including the value of km(y^) 
depends only on the geometry of the neighbourhood of the boundary S (see e.g. 
[41,42]). 

In summary, combining the two compatible bounds for edge states, we have the 
very constrained chain of inequalities 

- Km(T^b) Vh,«'h)m < ($S,-Am«>h)m <- KM(lfy) Vs ,$h)h • 

However, there is a significant difference between the two bounds. The upper bound 
shows the existence of edge states with negative energy. This is so because it means 
that there is at least one eigenvalue with negative energy. However, the lower bound 
applies to any state and establishes a global lower bound which guarantees stability 
of the vacuum for theories with mass larger than this bound. 

In the asymptotic regime p, > fj, ^ , the bounds become 

- (T^($s , ^h)m < (d>H , -Am4>h)m < - , ^h)h. (2.24) 

Since the inequalities hold for any value of 5 G (0,1) they become even more 
stringent, 

, ^s)m < ($H ,(2.25) 

in agreement with the asymptotic results of Refs. [41,42]. For small sizes y > ;^ 
p > fj,, the bounds become 

< ($s,-Am$s)m < -;^(d>3,$s)s . (2.26) 

This state becomes a zero mode in the case of Neumann boundary conditions (p = 
0). 
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The above inequalities show that negative energy states appear as long as p > 0 . 
The value of the energy gets progressively more negative the larger is the value of 
the Robin parameter, or equivalently the smaller is the manifold. The splitting 
procedure used for the proof shows that the lowest eigenstate can be approximated 
in the direction normal to the boundary by the lowest eigenstate of the 1 dimensional 
Laplacian with Robin boundary conditions. The latter, as it is easy to check, gets 
more localised at the boundary as the value of p increases. 

In summary, we have proved the following theorem. 

Theorem 2.2. Let M he an oriented Riemannian manifold with regular oriented 
boundary dM. The Laplace-Beltrami operator —A restricted to smooth functions 
satisfying the Robin boundary conditions (2.2) has at least one negative eigenvalue. 
As p = inf p becomes large the corresponding eigenfunction gets more localised on 
a small collar S around the boundary DM. 

The local character of the Robin boundary condition allowed us to find a lower 
bound in the spectrum. However, there are more singular boundary conditions 
which do not have such a bound for any size of the physical system [8,33]. 

3. Edge States of the electromagnetic field 

In pure electrodynamics on the space-time M x M the dynamical fields are one- 
forms A G A^{M) over the space M. The corresponding potential induced by the 
magnetic term of Maxwell action is up to a surface term 

where d* is the codifferential and (•, •) is the canonical scalar product of one-forms 

{a,l3)=[ g~^{a,l3)dpg, a, jS G A\M) , (3.2) 

J M 

induced by the Riemannian metric g of M. Here, g~^{- ,•) denotes the canonical 
scalar product of forms induced by g, dpg the Riemannian volume on M and 
* : A^(M) —)• the Hodge star operator, i.e. the unique operator that 

verifies 

g~^{a , (3)dpg = a A-kfd, Va,/3 € A^(Af) , (3.3) 

where n is the dimension of M. The following identities will be needed later on, 
cf. [48]. Let a,l3G A'=(M). 

{a,P)=[ g~'^{a,P)dpg= [ a A */3 , (3.4) 

JM Jm 

(a ,/3) = (*« ,*/3), *I = d^g, *dpg = 1, ★★ a = (—I)*^”“^^a . (3.5) 

In quantum electrodynamics the Hodge-de Rham Laplace operator 


-A = d*d 


(3.6) 
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plays the same role as the Laplace-Beltrami operator in scalar field theories. 

Using the Hodge star operator, the Laplace operator defined in Eq. (3.6) can be 
written, for A € A^{M ), as 


d*dH = (-l)"+i*d*dy4. (3.7) 

Gauge invariance implies that the Hodge-de Rham-Laplace operator (3.6) has 
an infinite number of zero-modes of the form A = dij). One can get rid of these 
spurious zero-modes by fixing the Coulomb gauge condition d*H = 0 or using 
Feynmann gauge by adding a term of the form dd* to the operator (3.7) to get the 
standard Hodge-de Rham Laplace operator 

-Ai =d*d-bdd*. (3.8) 

In the rest of this section we shall prove that for a particular choice of Robin-like 
boundary conditions the non-elliptic operator d*d is bounded from below indepen¬ 
dently of the gauge fixing condition chosen. 


3.1. Boundary term for Laplace operator on one-forms 

As in the previous section, we will use the quadratic forms associated to the oper¬ 
ators rather than the operators themselves to obtain the bounds. One important 
ingredient in the proof of Theorem 2.1 was the use of Stokes Theorem to identify 
the boundary term associated to the Laplace operator. We shall first derive Green’s 
identity for the Laplace operator (3.6) and write the boundary term in an explicit 
form. 

Using the identities above, one can get for A G K^{M ), 


(A,-AA) = (-1)”+W g-\A,i.d*dA)diLg (3.9a) 

Jm 

= (_l)"+i/ AA**d*dA= [ AAd*dA (3.9b) 

JM Jm 

= f dA A :*:dA — f d (A A *dA) (3.9c) 

Jm Jm 

= {dA,dA)- [ j*{AA*dA), (3.9d) 

JdM 


where j : dM ^ M is the canonical embedding of the boundary dM into the bulk 
manifold M and j* denotes the corresponding pullback map of forms. In (3.9c) we 
have used the identity 


d(A A *dA) = dA A *dA — A A d * dA 


and in (3.9d) the general Stokes’ Theorem, cf. [1, Theorem 8.2.8]. 
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Using standard algebra of imbedding maps [21], the boundary term in (3.9d) 
can be rewritten as 


[ j*{AA*dA)= [ j*AA-kQM3*\N<^A 

JdM JdM 


IdM 


9 ^{fA,fwdA)dpg 


g \a,a)dpg = {a,a)dM , 


(3.10) 


IdM 


where g is the pull-back of g to dM, -kQM denotes the Hodge star operator with 
respect to g, Af the vector field defined in a neighborhood of dM by the geodesic 
flow induced by the outward normal vector n of the boundary dM, and ij\f is the 
contraction with Af. Also we have defined a := j*A and a := j*\j 3 'dA maintaining 
the capital-small letter notation to distinguish objects on the manifold M and its 
boundary dM . Finally (•, ■)dM stands for the canonical scalar product among one- 
forms at the boundary: 

{a,l3)dM= / g~^{a,l3)dg,g = / aA-koM/3. (3-11) 

JdM JdM 

Compare the latter equation with Eq. (3.2). 


3.2. Robin-like boundary conditions and semiboundedness 

In spite of the fact that the Laplace-de Rham operator d*d is gauge invariant 
like the original potential V in (3.1), the boundary term (3.10) is not manifestly 
gauge invariant. Thus, the most general boundary conditions break gauge invari¬ 
ance. However, there are large families of boundary conditions which preserve gauge 
invariance, e.g. the Dirichlet boundary condition (a = 0,0^ = j*ijg'A = 0) or the 
Neumann boundary condition d = 0. One particularly interesting family preserving 
the gauge symmetry is given by Robin-like boundary conditions defined by 

rCf^iA - d(-Ao)-M*A) = Kj*{A - d(-Ao)-M*A) (3.12a) 

*dM * (A - d(-Ao)-M*A) = 0, (3.12b) 

where Cjg- is the Lie derivative with respect to the geodesic vector field N defined in 
a neighborhood of dM by the geodesic flow induced by the normal vectors n of 9M, 
Aq is the Laplace operator acting on scalars with Dirichlet boundary conditions and 

K :dM ^ ]g(n-l)x(n-l) 

is a function with values in the space of (n — 1) x (n — 1) real orthogonal matrices 
(as required by the reality of the electromagnetic potential A) with respect to the 
pointwise scalar product among one-forms defined by g ~^{-, •). 

The invariance of the boundary conditions (3.12) under gauge transformations 
<j) which vanish at the boundary j*(f) = 0 follows from the fact that in such a case 

d(j)-d{-Ao)-^d*d(j) = 0. (3.13) 
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Notice that since Aq is the Laplace operator acting on scalars with Dirichlet 
boundary conditions, the boundary conditions (3.12) reduce to 

fC^{A - d(-Ao)-MM)(p) = K{p) 3 *A{p), (3.14a) 

*3M *{A- d(-Ao)“M*7l)(p) = 0 (3.14b) 


for any p G dM. 

The most relevant property of these boundary conditions is the cancellation of 
the boundary term involved in the proof of self-adjointness of the Laplace-de Rham 
operator d*d. Indeed, such a boundary term can be rewritten as 


{A , -AA) - {A, -AA) = (a, d)dM - {a,d) 


dM 


/ j*AA*dMj*^N^— / j*A A *dMj*<iij^A 

'dM JdM 


IdM JdM 


(3.15) 


or as 


since 





fA A *9MfC^iA - d(-Ao)-M*7l) 
fA A *dMfdiMiA - d(-Ao)-MM) 
fA A *aMfCMiA - d(-Ao)-M*i) 
fA A *dMfdwiA - d(-Ao)-M*i), 


(3.16) 


a = fC^A - fdwA = fC^iA - d(-Ao)-M*7l) - fdi^{A - d(-Ao)-M*Al). 

Because of the boundary condition (3.14a), the contributions of the first and third 
terms in (3.15), reduce to 

(a, Kd)oM — {Ka , d)dM] (3-17) 

which cancels out due to the orthogonal character of K. The remaining terms 

[ fAA*9MfdwiA-d{-Ao)-f*A)-[ fAA*9MfdwiA-d{-Ao)-f*A) 

JdM JdM 

can be rewritten as 

- I {dl^fA)fiMiA-d{-Ao)-f*A) + f {dl^fA)fwiA-d{-Ao)-f*A) 

JdM JdM 

by integration by parts, where d*^ is the co-differential operator of dM. Now, 
since 


d*(A-d(-Ao)"MM) = 0, 


(3.18) 
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we have that 

j* * d * (A - d(-Ao)-MM) = - - d(-Ao)-MM) = 0, 

(3.19) 

because *dMj*A = — j*w * ^ 

j* * d * ^ = *dMj*^N' * A- dl^j*A. (3.20) 

The second term in (3.19) vanishes due to the boundary condition (3.14b). Then, 

d*J*A = 0, 

which implies that the two boundary terms (3.18) also vanish. In summary, the 
boundary condition (3.14) defines a symmetric extension of the Laplace-de Rham 
operator. It is easy to show that this extension is also essentially self-adjoint. 

We shall show now that the self-adjoint extension of the Laplace-de Rham op¬ 
erator (Eq. (3.6)) acting on one-forms constrained by these boundary conditions 
(3.14) is bounded from below. 

Theorem 3.1. Let M be an oriented Riemannian manifold with regular oriented 
boundary DM. Let n be the outgoing normal vector, N the veetor field defined in 
a neighborhood of dM by the geodesic flow induced by n and K a smooth function 
with values in (n — 1) x (n — 1) orthogonal matrices, defined at the boundary dM . 
The Laplace-de Rham operator — Am restricted to smooth one-forms A € A^(M) 
satisfying the Robin-like boundary conditions (3.14) is essentially self-adjoint and 
bounded from below, i.e. 

{A,-AMA)>-^il\\A\\f (3.21) 

where /Xq is a finite positive constant. 

Proof. The proof will be almost identical to the proof of Theorem 2.1 and will 
also be derived in several steps. Only Step ii) differs substantially. Again we shall 
use the splitting of the manifold into a collar neighborhood of the boundary 5 and 
its complement, cf. Fig. 1. 

i) The collar Laplacian. 

Let —Ah be the Laplace operator of Eq. (3.7) on the manifold S defined on 
the domain 

Vs = {Ag \j*CN{A - d(-Ao)-M*A)(p) = 0 ,p € , 

j*CAf{A - d(-Ao)“M*A)(p) = K{p)a{p),p G dM ; 
*dMj*C^ * (A - d(-Ao)-M*A) = 0,p G as} , (3.22) 

where dE = dM U dE_ and (S) is the Hilbert space of one-forms (cf. Eq. 
(3.4)) whose coefficients are in 77^(S). This operator is essentially self-adjoint, 
as shown in the previous subsection. 
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ii) Bound on the collar Laplacian. The collar Laplacian and therefore the 
associated quadratic form are bounded from below, i.e., 

where po is of the order of sup|itr(p)| and \K{p)\ stands for the norm of the 
matrix K(p) . Again the bar over the domain stands for the closure with respect 
to the corresponding norm. We postpone the proof of this step until the end of 
this section. 

ill) The inner Laplacian. 

Let —Am\s be the Laplace operator on the manifold M\S defined on the 
domain 

2?m\h = {a G nUM\E)\j*C^iA - d(-Ao)-id*A)(p) = 0, 

^aMj*C m* (A-di-Ao)-^d*A) = 0,pGdE_}. > 

This operator is positive and by [18, Theorem 7.2.1] its quadratic form has 
domain 'H\i (M\S) without any constraints at the boundary. 

iv) Equivalence of the bounds on the operator and on the associated 
quadratic form. 

By the same arguments used in Steps iv) and v) of the proof of Theorem 2.1, 
in order to find a bound for —Am with domain 

Vm = {AgHI. (M) IfCMiA -d(-Ao)-M*A) = KfA, 

*om3*Cm * (A-d(-Ao)-M*A) = 0} (3.24) 

it is enough to find the following bound for its quadratic form: 

(A , — AmA)m > —Moll^lP ) AgVm . 

For A e q 7{(^j(M\S) we have that 

(A, —Am A) M = (A, —AhA)h + (A, —Am\eA) 

> -^^l\\A\\l > -HlWAfM , (3.25) 

where we have used ii) and iii). □ 

Proof of Step ii). The proof will follow similar lines to those of Step ii) of 
the proof of the Theorem 2.1, but in this case we will bound the quadratic form 
(A, — AhA)h by that of a one-dimensional problem with a flat metric. 

First we remark that the operator — A^ has an infinite dimensional kernel. More 
concretely, all the exact one-forms of (3.22) are in the kernel of the Laplacian. 
That is, if A^ = dcf) G then —AA^ = 0. One particular class of such exact one- 
forms A^ = d(/) are those defined by fields (p vanishing at the boundary of M, i.e. 
JdM^ = 0. The last condition guarantees that the zero modes satisfy the boundary 
conditions defined by (3.22) since = 0 and (3.13) implies that the boundary 
conditions (3.14) become trivial identities. 
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Moreover, the gauge field defined by 

A' := A- d^ 

satisfies the same boundary conditions as A because of the identity (3.13). 

Because of this freedom we can chose (p in such a way that A! is in Coulomb gauge 

d*A' = 0. This can be achieved for any gauge field A € by defining 

p=-Ag^d*A, (3.26) 

where —Aq is the self-adjoint Laplacian operator with Dirichlet boundary condi¬ 
tions densely defined on >C^(S). 

The Coulomb gauge fixing condition induces a decomposition of the Hilbert 
space £^(A^(S)) of gauge fields into two orthogonal complements 

/:2(a1(s)) n {H, d*A = 0} © /:"(A1(S)) n{A = d<p;<pe nl, (s), = o}. ( 3 . 27 ) 

Now, for any A G we have 

-AhA = -AhA' =- j^^drVdetgdrA' + d^d^A', (3.28) 

y^detg 

where dg and d^ denote the exterior differential and its adjoint in the Riemannian 
submanifolds r = cte. of S, and where we used the decomposition of the Riemannian 
metric in terms of the Riemannian normal coordinates on the collar neighborhood: 

= (orHc^)) ■ 

Notice that in this case the geodesic vector field Af is given by the tangent vector 
to the radial geodesics, i.e. = dr and the boundary conditions on the fields A' 
become 


drA'{p) = 0 , forp G fc- , 

fdrA'ip) = K{p)a'(jp), dr^/^A'^ip) = 0,forp G dM. (3.29) 


Thus, we have the following inequalities, 


1 


-(A, AhA)h = -(A', AhA')h = -(A', ^^9,yd^a,A')H+ |ld,A'f 

yjdetg 

>{drA',drA')E-{a,Ka')oM>{drA',drA')s-K{a,a')QM, (3.30) 


where K = sup|Ar(p)|. Let A^ax (p) > Amin (p) be the highest and lowest eigenvalues 
of g~^{p), respectively. Moreover, since the collar manifold S is compact, there is 
an infimum and a supremum of Amax (p) > Amin (p) in S. The same argument holds 
for the induced metric at the boundary g{RQ,9). The supremum and infimum of 
the eigenvalues will be denoted in this case by Amax and Amin • We can now get the 
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bound: 

P n—1 « n—1 

-{A , Ah21)h > Amin / V' - ^Amax / V' ajOid^lg 

>Amin(l-5)/ £[/anA'5nAl'dr-7T^%f^|a'Pld/ig. (3.31) 

The term in brackets corresponds for each i = 1,... ,n — 1 to the one-dimensional 
problem analogous to the one appearing in (2.11) and thus it is bounded from below 
by the same constant as that of the operator of Eq.(2.12). Hence, 

n 

-{A , AhH)e: > -Amin ^ H ( A,5f(l-<5) ) W^iWs 

i=l 

^ (aSS)) P'lll = , (3.32) 

where 

C'='^(AflE?J))- (3-33) 

and we have used that ||^'||| > Amin ll^illl- Now, since ||a 1|P = ||A'|p -|- ||d())|p 
and C > 0 we also have that 

{A , -AhH)b >-C\\A III . (3.34) 

i.e. —C is a lower bound of the operator —A^ 

Due to the orthogonal decomposition (3.27) the lower bound (3.33) also holds 
for the restriction —A to the domain of gauge fields which satisfy the Coulomb 
gauge condition 

d*A' = 0. 

—A is also essentially self-adjoint in the domain of smooth gauge fields in the 
Coulomb gauge which satisfy the Robin boundary conditions (3.14) that now be¬ 
come 

fC^A = Kj*A , *9Mj*C^ *A = 0. (3.35) 

The only difference is a change in the degeneracy of eigenvalues and in the 
absence of trivial zero-modes. 

The analysis of the Hodge-de Rham operator —Ai of Feynmann gauge (3.8) 
requires an extra boundary condition of Robin type in A = Xd*A, where A is an 
arbitrary smooth function of the boundary dM. 

A similar analysis should lead to the existence of edge states and a lower bound 

{A ,—/S.\jA) >—^\\\A\\'^, (3.36) 

where the constant is also A-dependent. An alternative proof could be obtained 
by using the Weitzenbock formula and a lower bound of the space curvature. 
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4. Edge States for the Dirac Operator with Atiyah-Patodi-Singer 
boundary conditions 

The existence of edge states is not an exclusive property of scalar and vector fields. 
We shall show in this section that the Dirac operator with APS boundary conditions 
also has edge states. 

Before we start the proof for the existence of edge states in this case, we will need 
to introduce some definitions. We assume that M is a Riemannian spin manifold 
with a given spin structure. Let ir : E —>■ M he the vector bundle over M associated 
to the fundamental representation of Spin(n) and oj the spin connection of E. The 
covariant derivative of w in if defines a map 

V : r{E) X X(M) ^ r(E) , 

where r(if) stands for the space of sections of the vector bundle and X(M) is the 
space of vector fields on M. 

In Riemann normal coordinates, the Dirac operator acting on the restriction of 
the vector bundle tt : if —>■ M to the collar neighbourhood S = dM x {R — e, i?} of 
DM [12] reads 

0 =. Vat - Z 7 e • Ve + , (4.1) 

= -ijf\/j^ + A'{m) (4.2) 

where Jf = 7 ^A/)i with the components of the vector field J\f while —170 • Vg 
denotes the tangential component of ij) . We use the notation 

A'{m) = —ijeVe + 7717 '^^^ . (4.3) 

Let the restriction to the boundary of this operator be 

A{m) = A'{m)\r=R = -ije ■ Vg + 7717 '^+^ . (4.4) 

The matrix 7 *^+^ is hermitean, anticommutes with all the 7 ’s and (y^+i)^ = I. 
The Dirac operator above is purely spatial and does not involve any time derivative. 
A{m) is an essentially self-adjoint operator acting on the Hilbert space of square 
integrable sections induced at the boundary that we denote by 'H(i9M). This is so 
because we assume that the boundary dM is a smooth compact manifold without 
boundary. 

Let us define the Atiyah-Patodi-Singer self-adjoint extensions of the Dirac op¬ 
erator ]p. Performing integration by parts, we get Green’s identity for the Dirac 
operator, cf. [14], 

($,04-)-(0$,4') =i((p,#) 0 M , (4.5) 

where as usual we denote the restrictions to the boundary with small size Greek 
letters and ^ = 7 ^n^ with the components of the normal vector field to dM. 
The boundary conditions must force the R.H.S to vanish identically in order to 
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define a symmetric extension. Let JC be any self-adjoint operator on T-L{dM) such 
that 

> 0 (4.6) 

and 

ifiJC = -K.il. (4.7) 

Since K^ is strictly positive by condition (4.6), the Hilbert space at the boundary 
admits an orthogonal decomposition in terms of (dM), the positive and negative 
eigenspaces of K respectively. The second condition, Eq. (4.7), now guaranties that 
in+idM) = n-{dM). 

The generalised APS boundary conditions consist in requiring that ^\dM G 
T-L~ {dM) and therefore in defining the domain as 

S(|)) = {4- G H{E)\^9m G H-{dM)] , (4.8) 

where the Dirac operator is essentially self-adjoint, cf. [9,14]. 

In order to analyse the existence of edge states we will need to analyse the 
relation between ij) and another second order differential operator that one can 
dehne on the Riemannian vector bundle. The second covariant derivative is defined 
on a generic tensor field T as 

■= Vx'^yT — 'KxxyT . 

Definition 4.1. The connection Laplacian is the operator defined by 

-A := -V^V^ . 

There is an important relation between this Laplacian and known as 
Lichnerowicz-Weitzenbock identity, cf. [36]: 

# =-A-hR-hm^ . (4.9) 

where R is the Ricci scalar curvature of the Riemannian metric of M. 

This theorem defines a formal relation between the operators. Now we want to 
analyse the existence of edge states for the operator Ip and this will be equivalent 
to analysing the existence of edge states for the operator Ip . But we need to select 
a domain for the latter that is compatible with the domain of Ip. We will define 
this domain by 

©(02) = g D(.0)|04' G S(0)} . (4.10) 

With this choice, we have that the following identity holds for all 4*, $ G T>{p'^): 
(4>, 0^4-) = (0$ , 04-) = (02$ ^ . 

The boundary terms vanish by means of the conditions in (4.10). In fact, a general 
argument using the spectral theorem shows that D{lp ) is an essentially self-adjoint 
domain for Ip. 
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We will now make a particular choice for 1C and hence of the boundary condition. 
We will show that this particular choice leads to low lying edge states. 

As the operator K. for the generalised boundary conditions, we take the operator 

/C(/i) = . 

It is clear that it verifies conditions (4.6) and (4.7). If we denote by P+ and P_ 
the orthogonal projectors onto the boundary subspaces 'H'^(9M) and 'H~{dM) 
respectively, the boundary condition on the domain of 0 can be expressed as follows: 

S(0) := {^ll g H{M) I P+i; = 0} . (4.11) 

Let us write explicitly the boundary condition for JjP . We have that g 
implies that G ®(0) • 

Now 

+A'(m)4' (4.12) 

{$'^)\dM = -iiP + A{m)ip , (4.13) 

where ip = (VAt4^)|aM and ip = 'i’\dM ■ According to (4.10) and (4.11) the domain 
of Ip^ consists of the functions that fulfil P+ip = 0 and P+{lp’i’)\dM = 0, which 
according to the identities above can be written as 

P+ip = 0 

P+ —iipip + A{m)ip =0. (4-14) 

Recall that P^ is the projection onto the positive (or more generally non¬ 
negative) space 'H'*’(9M) of the operator JC{p). Our aim is to show that there 

exist low lying edge states. In the spirit of the proof on the existence of edge states 

for the Laplace operator given in [8], we are going to select an edge state satisfy¬ 
ing the boundary condition above. In order to do so, we need to analyse first the 
relation between the operators JC{p) and A{m). First notice that 

A{m) = A{ii) + (m — , (4-15) 

A(/i) = {ip)‘^A{^) = —ipLlC{iP) . (4.16) 

Hence the boundary condition (4.14) can be rewritten as 

P+ —ii/iip — ii/ilC{p)ip + {m — pPy'^'^^iP =0. (4-17) 

According to this, the most simple case occurs for m = p . Therefore, for simplicity, 
consider from now on that p = m and that the boundary condition is 

P+lC{m)ip = 0 

P+ —iifiip — ipilC{m)ip = 0 . 


(4.18) 
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Moreover, there is a direct relation between the square of the operator K,{m) and the 
connection Laplacian on the boundary of the manifold given by the Lichnerowicz- 
Weitzenbock identity, see Dehnition 4.1 and Eq. (4.9): 

/C(m)^ = i'jiA{'m)i'jiA{'m) = = A{m)'^ = —Aq + + Rom ■ (4.19) 

Now Rqm stands for the scalar curvature of the induced Riemannian connection at 
the boundary. 

In contrast to the situations on Sections 2 and 3, the operator is automati¬ 
cally positive and thus there is no need to prove its semiboundedness. But we will 
show that there are edge states, whose eigenvalues are close to zero and therefore 
below the mass gap, which is not obvious. We will follow the approach taken in Sec¬ 
tion 2.2 and find an upper bound for the energy of certain states that are strongly 
localised at the boundary. 

First we select an eigenstate ^ € R- of /C(m) with eigenvalue —A, with A > 0, 

/C(m)C = -AC . (4.20) 

For this state we can prove the following bound. 

Lemma 4.1. Let C he an eigenstate of Kirn) with eigenvalue —A and let —Aq 
be the boundary Laplacian of the spin vector bundle induced at the boundary dM. 
Then 

(C,-AeC)<(A 2 -m 2 )||Cf+ a„,ax||Cf , 

where CTmax := supp^g^ I^9 m(p)| • 

Proof. 

From (4.19) and condition (4.20), we have that 

(C 5 {-Ae + 171^ + TZdM)0dM = A^IIClP 

(C )-AeOdM = (A^ - m^)||C|P - (C > RomOom 
^ (C)-AeC)aM < (A^ - m^)||C||^ + cTinaxllCII^ ) (4-21) 

□ 

Notice that we can choose A such that A^ — -I- Umax is the closest element 

to the smallest eigenvalue of — A^. It is well known that the eigenvalues of the 
connection Laplacian, cf. Definition 4.1, on a manifold without boundary decrease 
with increasing volume of the manifold. As an example consider the situation for 
the sphere S'" with radius p. The eigenvalues of the Dirac operator are given by 
A; = ±(/ -I- nl2)lp , cf. [16]. Using the Lichnerowicz-Weitzenbdck identity and the 
fact that in this case the scalar curvature is i? = n{n — l)/(4p^) one gets that the 
lowest eigenvalue for the connection Laplacian is n/(4p^). Hence, by increasing the 
volume of the manifold it is possible to make the smallest eigenvalue different from 
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zero arbitrarily small. 


Now we will show (see Theorem 4.1) that if — rn? + dmax is small enough 
and m is large enough, for a given width e of the collar neighbourhood S, there 
will appear negative energy states for —A, which means that there will appear edge 
states for 0^ with energies close to zero. Since 

(0^',04') = (^',0^4') 


for 4' fulfilling (4.10) there are corresponding low-lying edge states. In practice 
these conditions will be met, as stated in the previous paragraph, for large enough 
manifolds. 

Consider the following vector state with ^ taken as in Lemma 4.1: 

4* = exp (^fctan(^(i?o - r))) ^(0) , (4.22) 

for i?o — e < r < i?o and vk = 0 if r < i?o — e- Since 4'|r=i?o = 'tp = ^ , it is clear that 
'k € ■ Now we need it to be in the domain of and therefore we need to 

select k in order that the boundary condition (4.18) is satisfied. It is enough that 


tp = —IC{m)'ip . (4.23) 

A short calculation shows that for tk as in Eq. (4.22), the left hand side is i/f = 
“IfC while, according to (4.20), the right hand side verifies —IC{m)tp = A^ . Hence 
4- e if 

k - 

TT 

Moreover, under the latter condition the state satisfies 


i) = kxP . (4.24) 

Let us now compute the mean value of p^ in this vector state. We can assume 
that on the collar neighbourhood, we work on a frame of the vector bundle E such 
that Vj^'k = dr'^. Then, 


P 

(4-, -A^-) = (V^-, V^) - A(^, OoM = / / 

J Rn — e J d 


fRo r 

9^' 

/ / d/Xg 

JRo-e JdM 

dr 


rRo 

>Ro-e JdM 

< (1 + '^)IICIIaM 

r-fl, 

iRo-e 


2 

dM 


[ [ d/igexp (2A:tan Alien 

JRn-eJdM \ Ze / 

(A^ -Rp + CTmax) f exp (2fctan ;^(i?o - r)) dr 

JRn-e \ Ze / 


r "^0 / TT \ 

-b / A^(I-b tan^( — (i?o — T")))^ exp ( 2fctan — (i?o — r) 1 dr 
JRo-e 2e V 2e / 

< (1 + '^)lieilaM A (^- -b + (A^ — -b crmax)e 


A 


l + d 


2 I6e2A2 1 -b S 
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Here we have used the explicit form of the function 'I', Eq. (4.22), Eq. (4.24), 
Lemma 4.1, the identity 

dt(l + tan^ t)^ exp(—2A:tant) = ^ H—(4-25) 

2k 

and 6 is defined in Eq. (2.10). Thus we have proved the following: 



Theorem 4.1. Let —A be the eonnection Laplacian, see Definition 4-1, defined 
on the Hermitean bundle tt : E —)■ M and let 'I' € S)(.0^) be a state of the form 
shown in Eq. (4.22), with f{9) satisfying Eq. (4.20) and where k = Then the 

following bound holds. 


(4/ 




A 


1 TT^ 

2 16e2A2 



+ (A^ - + crmax)e 


Hence, as stated above, A^ can be chosen of the order of — (Tmax when the 
mass is large enough and if the lowest eigenvalue of —Ag is small enough. In this 
case, the first term in the brackets of the last inequality above becomes negative, 
showing that — A has negative eigenvalues. Moreover, for large A they become edge 
localised as follows from the boundary condition ip = Aip . In general, states with 
higher angular momentum might have energies greater than mf even though they 
are edge localised. 

Bulk states can be defined as elements of TdE), i.e., sections with compact 
support in the interior of the base manifold. Because for these sections '^\dM = 0, 
one has that 


(4',-A4') > 0 , 

and therefore 

(4-,02^) > rn^W'^f + (4',i?4') . 

In contrast to the situation for scalar fields, there might be no eigenvalues of Ip 
close to zero if the smallest eigenvalue of the boundary Laplacian —Ag is not small 
enough. Since this value depends on the volume of the boundary manifold, there 
will be a threshold value for the volume, which can be determined using Weyl’s 
law, cf. [38], such that edge states always exist. 

The difficulty of finding a negative upper bound for the expectation value of 
—A in the state (4.22) suggests the existence of a critical size threshold for the 
existence of edge states. We shall show in the next section that this occurs in some 
cases. However, to get a general result we will need a general lower bound on —A 
with the required transition properties. 

On the other hand, from the phenomenology of macroscopic devices the exis¬ 
tence of such a threshold for edge states is expected. For instance, for small enough 
samples of topological insulators the effects of low-lying edge states are not present 
and a threshold for their appearance is detected experimentally [35]. 
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5. Edge states of Dirac operator with chiral bag boundary 
conditions 


In addition to the former, there are other families of boundary conditions that one 
can consider, for instance, the boundary conditions used in the analysis of quark 
confinement [17] or their generalizations like the chiral bag boundary conditions, 
cf. [45]. These can be introduced using a general approach that leads to self-adjoint 
Dirac Hamiltonians. 

If the boundary term (4.5) is written as the difference of the two chiral compo¬ 
nents dtj- = i(l ± 7 ' n)4' of spinors, -|- , then we get 

S(4',$) =i(4'+,4>+)(5.1) 

It is easy to check that boundary conditions of the form 

{1 — Y • n)'>P = U{1-\-Y • n)'4’, (5.2) 

where 17 is a unitary operator on the boundary Hilbert space of spinors commuting 
with 'y ■ fi , make the boundary term vanish. These boundary conditions define 
domains of essential self-adjointness provided that certain regularity conditions on 
the unitary operator U are satisfied, cf. [8,15,32]. 

We will assume “local” boundary conditions imposing that U is a, finite dimen¬ 
sional matrix acting only on spinor indices. Using the Cayley transform, we can 
express (5.2) as 

1 _ [7^+1 

j + C/y+i V^- (5-3) 

As a particular case, we can consider U of the form 

[/ = e2*a''ctane“^ (5.4) 


which because of the identities 
l-Uyd+l + U 
I + Uyd+^~I ~ 


d-l-l 


1 - u 

I + U 


(l + 7‘^+i) 


=f cot(arctane“)(l — 7'^'*’^) — i tan(arctane“)(l -I- 7“^^^) 
='ie"“(l - 7^^+^) - ie“(l -b 7^^+^) = , 


corresponds to the chiral bag boundary conditions [45,13]: 

i ^1 — 7 • f/' = 0 • 


(5.5) 


Example: The 3D-Ball B 3 

Let us consider a Dirac electron moving on a 3D ball B 3 of radius Rq. The Dirac 
Hamiltonian 


H = —iy ■ V + 7ow 


(5.6) 
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is subject to the boundary condition (5.5): 


(5.7) 


with 7 • n = 7 i sin 9 cos ip + D 2 sin 9 sin (/? + 73 cos 9. 

The Hamiltonian H is essentially self-adjoint in the space of smooth functions 
satisfying (5.5). Let us consider the spectrum of H. It is invariant under rotations 
with generators given by the total angular momentum J* = L* -|- S'* where S* = 
— |eb^{ 7 j, 7 fc} (i=l,2,3). Let us consider stationary states of the form (see e.g. 
[24,46,47]) 


where, for k > 0 , 


'^37n{r,9,p) = 


^j,m,K{9,'p)(t)i{r) 

i^j,m-K.{9,p)(l>2{r) 


1 / Vj + IP 


and 




V j - TO + 1 y ™ i 


V 2 J +2 


(5.8) 


are the bispinor eigenfunctions of the three commuting operators = 

7 o(2S*Li -I- 1) associated to angular momenta. 


jm{r,9,p) = j{j + l)'^j^{r,9,ip) 

J3'^jmir,9,ip) = m'^jm{r,9,p) 

K^jmir,9,ip) = K'lijmir,9,ip) = (j + 5 ) 'i'jmir,9,ip). 


1 


Now, since 

-ij ■ V'^jm{r,9,p) = -ij -n + ^(1 - 7 o^^) ) ^im(r, 6 », 7 )), 

the solutions of the Dirac equation 

H^jmir,9,ip) = E^jm{r,9,p) 
must satisfy the coupled equations: 

dr + ^-{E + m)(j)2{r) = 0, 


(5.9) 


(5.10) 


— dr — 


1 -f K 


Mr) + {-E + m)())i(r) = 0, 


(5.11) 
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which can be decoupled into a pair of second order differential equations: 


di + -dr - 

r r 


(k- 1) 


dl + -dr - 


k{k + 1) 


r 

Since {7 • n, if} = 0 and 

7 • (/j) = 


(jjiir) = (m^ - 

hir) = [m^ - E'^)(l32ir). 

-nj^rn-K{d,(p)4'l{r)) ’ 


the boundary conditions (5.7) imply that 

<('2(i?o) = e“(()i(i?o). 


(5.12) 

(5.13) 

(5.14) 

(5.15) 


The solutions of (5.12) and (5.14) with smooth behaviour at r = 0 are given in 
terms of the Bessel functions [2,40], 


1.0 

0.5 

0.0 


-0.5 



X 

Fig. 3. Charge density distribution ilitvi/ of edge state with zero energy on a 3D ball B 3 for 
a = 1.16144 and mRo = 1. 


y/m + E _ 


y/m — E 


Mr) = "V "" 4 -} , (j) 2 {r) = — I. 


«+i , 


(5.16) 
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where the normalization was chosen to match the Dirac equations (5.10) and (5.11). 
Now, the boundary conditions (5.7) imply that 


(i?oVm2 - E^) 


{Rgy/m^ - E^) 


which is the spectral condition. The edge states correspond to states with negative 
kinetic energies (T = \E\ — to), i.e. —to < E < m and the simplest one corresponds 
to j = 1/2 and k = 1 , i.e. 



These states always exist provided that a > 0. 

Notice that in the limit a -A oo, (5.15) leads to Dirichlet boundary conditions 
for (/i, i.e. ())i(i?o) = 0 whereas in the limit a -A —oo we obtain Dirichlet boundary 
condition for (j) 2 , i.e 4 ' 2 {Ro) = 0 and because of (5.11), Neumann boundary condition 
for (j)i when j = ^, i.e. (j)[{Ro) = 0. 

In particular, we have a zero mode for a = log Is (Rom) — log 1 1 (Rom) which 
corresponds to the maximally localised edge state. The concentration of the state 
on the edge increases as the mass gap increases, which provides the perfect situation 
for a topological insulator. 

One can find more states with higher angular momenta and negative kinetic 
energies. And whenever a = log I(Rom) — log Ij_i {Rom) , we have zero modes. 

Notice that because the boundary condition (5.7) preserves rotational invari¬ 
ance, the energy spectrum is degenerate. For angular momentum j , there are 2j + 1 
degenerate states. There is an additional symmetry of Dirac equation in If we 
interchange the two bi-spinors in the ansatz 




/Dj,™ _„(6», (p)(/i(r)\ 


(5.19) 


it is possible to find solutions with the same energy. However, in the ball the 
boundary condition breaks this symmetry, and in particular there are no edge states 
of the form (5.19) when a < 0. This is in agreement with the experimental result 
that for small enough samples of topological insulators the effects of low-lying edge 
states are not present [35]. 

The number of such states is always finite and depends on the mass gap. For 
larger masses, there is a larger number of edge states. For mass to = ^, there 
are no edge states for boundary conditions with a < 0.4. The number of energy 
levels leading to edge states for different boundary conditions a = 1,2,3,4, 5 is 
3,11,18,35,98, respectively. 

A final remark concerns the breaking of chiral symmetry induced by edge states. 
Since chiral boundary conditions break this symmetry, there is an induced asymme¬ 
try in the spectrum between positive and negative energy modes. In this sense the 
particle-antiparticle symmetry is broken by the effect of chiral boundary conditions. 
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Fig. 4. Number edge states for different boundary conditions a. = 2,3,4, 5 on a 3D ball Bz for 
m/?o = 1- 


The above spinorial edge states do have a special meaning in QCD in a three- 
dimensional ball with chiral boundary conditions. The MIT bag model uses the 
chiral bag boundary conditions in the limit a —> oo. In that case, there is an 
infinity of edge states. In particular the lowest energy state is an edge state. The 
states of pions and protons made of quarks localised at the edges of the bag is a very 
interesting picture of the nucleon where according to asymptotic freedom, quarks 
will move quite freely inside a hadron. 
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